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Abstract 



Nowadays it is widely accepted tiiat the evolution of the universe was driven by some 
scalar degrees of freedom both on its early stage and at present. The corresponding cosmo- 
\ logical models often involve some scalar fields introduced ad hoc. In this paper we cultivate 

■ a different approach, which is based on a derivation of new scalar degrees of freedom from 

fundamental modifications of Einstein's gravity. In elaboration of our previous work, we here 
investigate properties of the dilaton-scalar gravity obtained by dimensional reductions of a re- 
cently proposed affine generalized gravity theory. We show that these models possess the same 
Oh< symmetries as related models of GR with ordinary scalar fields. 

As a result, for a rather general class of dilaton-scalar gravity models we construct additional 
first integrals and formulate an integral equation well suited for solving by iterations. 



(N 

^ ■ 1 Introduction 

Q>^ ' The cosmological observations of recent decades revealed that the universe expanded with accel- 
■ eration on two different stages of its evolution: in the very beginning and at present time. In 
CN . addition, the presence of new sorts of matter, the amount of which significantly exceeds that of 
\ the known matter, is also established very well. The models within the current paradigm of the 
Priedmann-Lemaitre- Robertson- Walker cosmology and the Standard model physics can hardly de- 
scribe all of them. This probably means that partial improvements of the existing theory may 
prove to be insufficient and the entire paradigm should be changed [IJ. 

For example, one can try to extend the particle physics by adding inflaton or quint-essence 
. fields, or consider the modified gravity, or use the paradigm of string/brane world (see, e.g., 
[21 El m [5] as a drop in the ocean of papers and reviews on the topic). 

Incorporating new scalar fields into the Standard model could be a technically simplest solution, 
but this approach lacks a fundamental basis. Moreover, it can hardly be falsified, since any 
disagreement with the observational data can be cured by adding more fields to better suit the 
data. Actually, now it is widely accepted that the string theory, being a fundamental one, suffers 
the same problem, with numerous compactifications which can provide any effective theory. 

Among the modifications of gravity there is a very interesting theory, originally introduced by 
Weyl, Eddington and Einstein, which has a fundamental distinction from General Relativity. It 
establishes the priority of geodesies within non-Riemannian geometry, while the metric arises at 
the level of the effective theory. On the other hand, the constraints on the initial theory come 
from the requirement to reproduce the metric theory of GR gravity. The constraints can be 
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resolved, and the use of a special two-step Lagrangian formalism allowed Einstein to formulate 
the effective theory with normal curvature term, cosmological constant and some 'dark' matter 
fields interacting only with gravity. These fields looked quite abnormal, and since the actual goal 
at that time (around 1923) was to unify the gravity and electromagnetism, the theory was not 
investigated in detail, soon abandoned and completely forgotten. 

Nowadays, the presence of dark fields is the fact of our life, and such a natural, almost in- 
evitable, presence of the dynamical cosmological constant is welcomed very much. Indeed, one 
of the intriguing questions of the dark energy problem is: why does it look so much like A-term 
[3 El El E] I while the origin of a so small A-term is very unlikely in the framework of GR and Stan- 
dard model. Thus it looks quite natural to re-establish investigations of the affine generalization 
of GR |im which can be called the 'affine' gravity. 

One of the basic features of the affine gravity is the presence of a massive (normal or tachyonic) 
vector field. In cosmology, the scalar fields are usually considered, but recently [11] we revealed 
the duality between this vector field (vecton) and the corresponding scalar field (scalaron) in case 
of reductions to two-dimensional theory. This allows us to investigate the basic properties of 
practically interesting cosmological and static solutions using the same approach as for gravity 
theories with scalar fields. In fact the affine gravity produces a scalaron field which significantly 
diff'ers from the standard scalar field in GR, thus compelling us to use a significantly generalized 
description of the dilaton-scalar gravity (DSG). And the obtained results can be applied even 
outside the context of affine gravity. For instance, the dimensional reductions of Yang-Mills fields 
[121 [ISl O [IS] also provide rather special scalar degrees of freedom. 

On the other hand, DGS theories constructed for the cosmological purposes, are usually rather 
simple. In case of affine gravity theory, the corresponding effective DGS model, that we call the 
dilaton-scalaron gravity (DSG), is defined rather strictly. Actually it looks quite cumbersome, 
being significantly non-linear in dimensions D > 3 and non-integrable even in the simplest linear 
case. Nevertheless, in this paper we will investigate the basic properties of the DGS model and 
try to answer the question: is the DSG model really complicated? 

The paper is organized in the following way. First, we discuss models which possess some 
fruitful properties, such as the presence of linear symmetry, the availability of iterative integration 
procedure, the existence of explicit non-trivial solutions. So we specify the set of simplest models 
within the class of DGS, mostly related to dimensional reductions. Then we show that DSG 
belongs to this set, despite its specific non-linearity. We thus conclude that affine gravity theory, 
being quite awkward at first sight, shares with common scalar field models in GR some useful 
properties. In Appendix we present some extra results which were derived while investigating 
the DGS gravity but can be applied to a rather general class of dynamical systems appearing in 
theories with vanishing Hamiltonian constraint. 

2 Dilaton-gravity theory with scalar fields 

A fairly general higher-dimensional field theories can be reduced to two-dimensional effective DGS 
models by taking into account their space-time symmetries. We specify the following generic action 

4ff =^/^[v'«(5) + W('p)(V^V')2 + Zi,((/p;V')V^V'iV%+X((^;V')] , /i = (t,r), (1) 

where (p is the dilaton field, and ipi denote all scalar fields of any origin, including the geometrical 
one. Note that in dimension D = 2 all the fields are dynamically equivalent to scalar fields or, in 
case of the massless vector, to the corresponding field charges. A practical transformation of the 
vecton into the scalaron in the dimensionally reduced affine gravity model was given in |11J. 

The most concrete results for these models can be obtained after the further reduction to the 
one-dimensional dynamical systems, describing static and cosmological states. For this, the system 
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is parameterized by a single coordinate, r, which can be time-hke or space-hke. The signs of the 
one-dimensional Lagrangian components depend on the parametrization, but they usually may be 
consumed by the free parameters of the model. 

So, the corresponding one-dimensional Lagrangian can be written in the following rather gen- 
eral form [m da da HBl [19]: 



(2) 



Here /i is a metric function, chosen so that /i > is a cosmological state, and /i < is a static state. 
The presence of the remaining metric gauge degree of freedom s, which is in fact the Lagrangian 
multiplier, leads to the vanishing Hamiltonian (constraint): 



H 



(3) 



After imposing this constraint we may fix s. For example, s = 1 is the lightcone gauge. In what 
follows we often replace the metric h with the new variable F defined hy F = h~^h + W{ip)ip . 
This reduces the dilaton-metric kinetic term to F(p and incorporates the Weyl transformation; the 
potential X(ip;tp) accordingly acquires the factor = exp(— j W{ip)dip). 



2.1 Linear symmetries and bilinear integrals in DGS 

Consider now the system ([2]) with just one scalar field -0 and the corresponding function, Z{ip), 
depending only on the dilaton variable. In this case, we may use a convenient gauge s = 1/Z{ip) 
and the new dilatonic variable ^ defined by = ip/Z{ip). Then our DGS dynamical system with 
the Hamiltonian constraint looks rather compact: 

C = Fi + i;'^ -e^U{^,ip), H = Pi + ijj^ + e^U{C,^) = 0, where U = nZX. (4) 

The corresponding Lagrangian system, 

i=-e^U, F = -e^U^, i^ = -e^U^/2, (5) 

actually contains only two independent equations, while the remaining one can be derived from 
the above Hamiltonian constraint H = 0. Yet it is more convenient to use all the three equations 
when searching for linear symmetries of the system. 

With this aim, let us introduce the coordinate vector = (F, ^, ip) and represent our dynamics 
in a more general form, 

C = Ajq'q^ -U{q), H = A.^q'q^ +U{q) = 0, (6) 
where in model ([1]): Ap^ = A^p = 1/2, A^^ = 1, and U{q) = e^U{£,,ip). The equations are 

2Aijq^ = -diU{q) ^ = (7) 

The last system can be contracted with a linear combination of the coordinates and some constants, 
{XAkj + Bkj)q^ + Cfc, where B^j is an arbitrary antisymmetric constant matrix, Ck is a constant 
vector, and A is a constant. Then we have: 

[{XAkj + Bkj)q^ + Ck\it ={[(A^fcj + Bkj)q^ + Ck]q''y - XAkjq^q^ = 

= --[{XAkj + Bkj)q^ + Ck\{A~^f'diA{q). 
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Using here constraint Q, H = 0, we replace the kinetic term, Akjcf q^, by the potential, —U{q), 
and conclude that the bilinear form. 



P = [{XAkj + Bkj)q^ + Ck]q' 



(9) 



is the integral of motion if the potential satisfies the linear partial differential equation: 



m{q) + -[Xq^ + {Bkjq' + Ck){A-^f']diU{q) = 0. 



(10) 



This theorem is applicable to the general constrained dynamics (gauge system) described by 
Eqs.([6]) with any constant matrix Aij and any function Z^(g^, (7")0 For our DGS system with one 
scalar field, the only nonzero components of the 3x3 inverse matrix A~^, are (A~^)^^ = (j4~^)^^ = 
2, {A~^Y^ = 1 and the potential has the simple dependence on the metric, U{q) = e^U{^,ip). 
In addition, in the three-dimensional coordinate space, the antisymmetric matrix Bi^j can be 
parameterized by a vector, Bkj = ekjihK It follows that Eq. ffTO]) splits into two independent 
equations for U{^,il)) (the second one expresses the condition of F-independence): 



This equation depends on seven parameters, but U is defined up to an arbitrary multiplier. Integral 
([9]), also defined up to a multiplier, depends on the same parameters and now reads: 



P = [(A/2 + h-i)i - h2il} + ci]F+ [(A/2 - h^)F + bitlj + C2] e + [Xtp + 62F - 61^ + C3] -0 . (12) 



Note that when &2 7^ and (A — 2^3) 7^ the potential defined by equations (jlip must be 
exponential in ip, i.e., U = m(i^) exp[(253 — X)tlj/b2]- Inserting this into the first equation we 
find that 'u(^) is also exponential (or constant). So the model reduces to the well studied class 
of multi-exponential DGS theories with the completely integrable subclass of Toda - Liouville 
theories (see [20] and references there) and we therefore expect to find in the exponential case 
some integrable systems. Below, we mostly suppose that 62 = (A — 263) = that does not mean 
ignoring exponential potentials. In fact, supposing that A = 6j = 0, it is not difficult to find the 
general solution of pT]) . 



where / is an arbitrary function and C4 — a new arbitrary constant. Making the simplest choice, 
f{x) = C5X we find that Y = giS^ + g2tp, with gi = C3C5 + C4 and g2C3 = cigi + 02- The potential 
depends only on the parameters 51, 32 and thus we are free to vary all the parameters q while not 
changing Y . For example, we can take c\ = 0, 2c2 = — (72C3 or C2 = 0, 2gici = —g2C^ and obtain 
two independent integrals discussed below. 

These examples show that relations between the potential and integrals are in general rather 
complex. We hope to demonstrate that the construction given by Eqs.Q - (jl2p is a powerful 
device, but it is not easy to use. In this paper, we only briefly discuss mathematics behind it and 
concentrate on simple but nontrivial examples. A very close but simpler approach to searching 
for integrals in DGS models can be found in |21j . 

^ Note that throughout the paper we consider the simplest DGS with one scalar matter mode. The general 
system is briefly discussed in Appendix A and multi-scalar DGS models - in Appendix B. 





y(e, i;) = \n\U\ = f{cs^ - 2ciV^) + CiC 



2C3 ^(C1C4 + C2) 



(13) 
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2.2 Master Integral Equation in DGS 



One may reduce the order of the system ([5]) by choosing ^ as the independent variable with 
dr = d^/xiO- The system of equations can be therefore simphfied if we define = F — Inx, 
Q = ip' , and use the combination of constraint Q with the first equation of the system instead of 
the second one. Then the system wih contain only the first order equations: 

^ = -e^U, T' = -g\ V' = a, {xQ)' = -e^U^/2. (14) 

Obviously, the solution for -0 and F can be written as formal integrals for Q and Q^: 

_ _ 

G{i)di^x{g-i], ^(e) = J-Q- / 

Now use these integrals in the remaining equations for x ^-iid xG and integrate them as well: 

'Co 



xiO -xo = - [ exp (X{g2;e}) UiC,l{G;md^^MG;C}, (16) 
2[x{0G{0 - XoGo] = - exp (X{g2;e}) U^{^I{g;^})d^ ^ Z2{Q\i]. (17) 



In the last equation replace x with the functional Xi and find the relation between Q and the 
functionals depending on Q and the potential U: 

r(c\ Xogp +X2{g;g}/2 , . 

We call this the 'Master Integral Equation', or MIE. It can be used to obtain the approximate 
solutions for the general potential. With a clever choice of the trial function for Q we can obtain 
a good approximation even with a few iterations. Although the integrations are in general rather 
cumbersome, the computations with MIE may be very effective near singular points or in vicinity 
of a horizon, where the solution can be expressed with the aid of convergent power series expansions 
(after explicitly resolving possible singularities). 



3 Partially integrable DGS models 

The models of gravity with scalar fields often can not be solved explicitly. Yet in case of just one 
massive scalar mode one may expect that deriving an additional first integral and/or applying 
iterative approximation procedure will allow to investigate the most important properties of the 
solutions. For example, in cosmology the slow-roll approximation technique provides a good infla- 
tion model easily derived for many non-integrable configurations. Having this in mind, consider 
some types of DGS configurations possessing the inner symmetries, which simplify their analysis. 



3.1 Integrable DGS models 

Let us first explore the completely integrable configurations. They may arise in case of the ex- 
ponential potential linear in its arguments: U = exp (gi^^ + §2^) ■ It possesses two commuting 
linear first integrals (the conditions for existence of several commuting first integrals can be found 
in Appendix). Indeed, such potential suits the condition (jlip with the non- vanishing parameters 
C2 = —cigi, or C2 = —035(2/2. The two independent integrals, 

Pi=F- gii P2 = 2iP- g2i (19) 
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are in involution (commute). Taking into account the constraint, we find that the number of the 
integrals is equal to the number of variables and thus the system is integrable. This configuration 
with linear exponential potential may be related to the Toda-Liouville systems [19j, or it can be 
derived from a quite special cylindrical reduction, as we will show below. 

Another example comes from the effective D = 2 massless scalar-vector configuration 



(^i?(g(2)) + W{ip){V^ip)^ + Z(v9)V^V'V^V + Y{ip)B^,B>''' + X{ip)\ , (20) 



where B^i, = d^A^ — d^A^. In the one-dimensional cosmological or static case, the scalar and 
vector fields can be integrated out, leaving as a trace the effective charge of the vector field, Qq, 
and the integration constant of the scalar field, cq = Z{^)^. So the dynamics may be described by 
the effective Lagrangian supplied with the Hamiltonian constraint: 

C = -Fi + e^U{i) + V{i), H = Fi + e^UiO + V{^) = 0. (21) 

Here the original kinetic scalar field term contribute into the effective potential V, while the kinetic 
vector term contribute into U. 

To derive an explicit solution of such a model we now need just one additional integral. Let us 
use here a sort of an inverse approach: first choose the suitable model which does possess such an 
integral, and then find a corresponding physical theory described by the model. First check the 
general condition (jlOp for the (bi)linear integral to exist. For the system ()2ip it can be satisfied 
only when one of the effective potentials, either U or V, vanishes, because the total potential 
term, e^U{^) + V{^), does not possess the linear symmetry. Therefore, the mixed system with 
scalar and vector fields cannot have such kind of integral. With both potentials non-vanishing, 
there remains a possibility for a non-linear integral to exist. There is no general algorithm for 
constructing non-linear integrals but, nevertheless, an integral quadratic in momenta was found 
for a dynamical system of this sort (in |16] . [17j). 

Reproducing the procedure of |16j in our notation, first rewrite the equation for the metric: 

F = -e^U'{0-V'iO- (22) 
For linear potentials, U{^) = gi^ + 53, V{^) = + 94, it becomes trivial and can be integrated: 

i^ = -5ie^-52, = C - gie^ - g2F = R{F). (23) 

Here C is the new first integral, quadratic in momenta and non-linear in the coordinate variable 
F. Multiplying it by ^ and expressing by using the Hamiltonian constraint we find: 

[e^(5iC + 93) + 92C + 94\F + 2{C - g^e^ - 52^)^ = 0. (24) 

This expression is linear in momenta and may be more convenient to deal with. 

Supposing that gi 7^ 0, let us rescale ^, ^ = 9i£, + 53, and rewrite the above equation as 

[{9ie^ + 92)i + 9m - 9293]F + 2i?(F)| = 0. (25) 
Since gie^ + 52 = we find the equation 

2RiF)i - iii + (3451 - 9293)F = 0, (26) 
which can be integrated after dividing it by R^^'^: 



2R~^/\F)l - 2i?-i/2(Fo)eo = (5253 - 9m) / R-^'\F)dF. (27) 
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It is instructive to rewrite the the final expression for ^{F) in the form: 

[r{f) 



RiFo) 



Co + ^{9293 - 9a9i)VR{Fo) [ R-^/\F)dF 



(28) 



With = 0, the integral in r.h.s. will contain arctan, arctanh or exponential functions depending 
on the parameters, otherwise it cannot be expressed in terms of elementary functions. Let us also 
note that the model with 92 = may describe the spherical dimensional reduction of the massless 
scalar and vector configuration [16^ I17j . For example, the Reissner-Nordstrom solution in the 
absence of the scalar field (when cq = 0) can be derived from (j'28|) when 92 = 94 = 0. 

We would like to mention that the full integrable model with all gi non-vanishing provides a 
rather rich dynamics as compared to the Schwarzschild or RN solutions. The investigation of such 
fairly complex integrable models can significantly improve our understanding of complex dynamics 
in non-integrable theories with massive fields. 

3.2 DGS and cylindrical reduction 

The general cylindrical reduction is, in fact, much more complicated than a spherical one. Consider 
the decomposition of the line interval in D = 4: 

dsl = {gij + if amn fT'P'j) dx'dx^ + 2ipi„^ dx'dy"" + ip amn dy"'dy'' , (29) 

where i,j = 0,1, m,n = 2,3, all the metric coefficients depend only on the x-coordinates {t,r)^ 
and y"* = (0, -z) are coordinates on the two-dimensional cylinder (torus). Note that ip plays the 
role of a dilaton and Umn (detain = 1) is the so-called (T-field. 

As was shown in [18], the reduction of the Einstein part of the four-dimensional Lagrangian, 
y—g^ R{g^^^), can be written as: 

,2 



(')) + ^(V^)' - |tr(VacT-iVaa-i) - 



(30) 



where (/?™ = difj^ — djipf. This complicated action can be simplified for some interesting cases. 
With a special choice of the matrix amn, 

a22=e^, 0-33 = 6"'^, O"23 = O"32 = 0, (31) 

and the vanishing y?^, one can derive the simpler effective action: 



(32) 



where Qi is the contribution of the integrated out Abelian gauge field 99^. This gives us a serious 
motivation to study the DGS models in which scalar potentials are exponential in ■0 . Their 
one-dimensional configuration can be described by a dynamical system Q, with the effective 
potential U = e52V'$(^). Actually, in pure geometrical case one has (p = e since ip/Z = ^ and 
Z{lp) = —p/2. The potential is power-like in p (even taking into account the Weil transformation), 
and thus exponential in ^: = e^^^. This is the integrable configuration considered above, but 
let us first consider a bit more general potentials with arbitrary ^(S,). 

It may be of interest to find how does the general function $(^) break the integrability. The first 
integral is obviously P = lip — g2C (or, equivalently, P = {g2 — 2Q)x using £^ as the independent 
variable). If P vanishes, we have the exact expression Q = 52/2 and then applying the MIE 
procedure easily derive the following special solution: 



V' = ?(e-eo), -^(0=-^o-f (e-Co), X(0 = Xo - e(^o-^'«°/^) f^me'"^"dt (33) 
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With P 7^ we return to Eqs. ()14p and show that the equation for x can be transformed into 
a simple Hamiltonian system with 'dissipation' vanishing of which transforms it into completely 
integrable system. Differentiating x' = —e^'^^'^^ ^{5^) and using the other equation and integral 

^' = a = (52-P/x)/2, J^' = -G^ = -{g2-P/x?/^, (34) 
we derive a somewhat unusual equation for x containing a 'dissipative' term: 

(In x')' = J" + 92i'' + (In c^.)' = ^ - + (In $)' . (35) 

This equation can be rewritten as a Hamiltonian system. If we denote x' = (^ = =^1), we 
find its canonical formulation, 

^''•"'^W "-f-i^^j +(>n*)' = -^, (36) 
with the ^-dependent ('time' dependent) Hamiltonian, 

H = H^^eeP xqII^ - x(ln (37) 

The total derivative of this Hamiltonian on the equations of motion is simply 

^ = -x(O(ln^(e))"/0. (38) 

For the integrable potential $ = e^^^, the Hamiltonian is conserved on the equations of motion, 
^ = 0, and one can derive the explicit exact solutions. In view of its simplicity, both technical 
and conceptual, this system may serve as a good laboratory for qualitative and quantitative studies 
of non-integrable DGS systems. 

3.3 DGS and spherical reduction 

The spherical dimensional reduction implies the following decomposition of the line element: 

ds\) = ds\ + ds\)_2 — 9ij dx^ dx^ + f'^'^ > (39) 

where v = {D — 2)^^ . It contains only the dilaton, so we consider the additional ordinary scalar 
field V' in the standard spherically reduced Einstein gravity: 



(40) 



The curvature term, ky rarely supports the symmetry of the scalar field potential ^{ip). Since 
we are investigating the contribution of the scalar mode, we omit the curvature term. Or, one 
may consider the vanishing kj^ as the result of the trivial cylindrical reduction. The resulting 
dynamical system is quite similar to the discussed above for the non-trivial cylindrical reduction, 
but it represents the other variant of the half-exponential potential: now U = e^^^^ilj). Indeed, 
Z{ip) = so ip/Z = ^ provides 99 = e^. Taking account of the Weyl transformation and of the 
gauge choice we have gi = u + 1, where normally = 1, 1/2, ... ,0, but we suppose that in this 
DGS model he parameter gi may be arbitrary. 
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In this model, there also exists a linear integral, P = F — giS^ = x{^' ~ 9i) + x'j but with 
the arbitrary potential this is insufficient for integrability. To find a partial result, let us 

combine first and last equations (|14p : 

If P = 0, we can obtain a simple first-order differential equation for y(^) = l/^(^). Recalling 
that x' /x = + 91 = G'^ + 91 and = dip/Q , we have: 

= + • (42) 

A nice explicit solution can be obtained if we take gi = 0: y = e~'^ / e"' dip. 

With non-vanishing P, simple equations with explicit Hamiltonian interpretation can hardly 
be derived. For instance, one can express Q = G{XtX') from the expression for P and put it 
into Eq. ()4ip . Then, again using the relation %[)' = Q one obtains a rather cumbersome second- 
order differential equation for x(V') with arbitrary ^{ip). Thus the system with the dilaton half- 
exponential potential looks more complicated than that with the scalar half-exponential potential. 



3.4 DSG and spherical reduction 

Finally we focus on the DSG model - a specific model of the DGS class, which is met in affine 
generalizations of gravity. The simplest case of DSG implies just one scalar mode, a scalaron tp, 
with the Lagrangian of the form [11]: 



if if 



(43) 



Compared to the ordinary scalar field, it has an abnormal coupling of the kinetic term to 
gravity. Now the dependence of ^-variable on the dilaton \s ip = since Z = \/ip instead of the 
standard Z ^p. In the gauge s = 1/Z , after the Weyl transformation with the factor O = p>'^~^ , 
we then obtain the effective power- like potential U{i,i}) = = ^'^/2-i$(V^7^). The function 

^("0^/0 is actually specified in the DSG model, but it depends on the space-time dimension, and 
for > 5 it cannot be given explicitly [10]. So it is much more convenient to treat it as an 
arbitrary function. For the ordinary scalar field potential, the dilaton field usually enters just as 
a factor. For the DSG, in contrast, the mixing of the dilaton and scalaron is essential. 

Again, we omit the curvature term in the potential since it does not correspond to the inner 
symmetry of the DSG model. The pure DSG potential, C/, satisfies the condition of the linear 
symmetry (|11|) . even in case of arbitrary $('0^/^), if we take A = 263, C2 = —I'h^ with vanishing 
all remaining parameters. Then the bilinear integral ()12p reads as 



P = iF-'^i + iPij. (44) 

Unfortunately, to find the other integrals or explicit solutions we have to require too much: P 
and 1^ must vanish and the function $ has to be specified. Instead, we may use another approach, 
the MIE, introduced in (jlSp . It suits the DSG system very well and allows to derive a power-like 
solution, as we will see now. 

We first choose Q = Cq^^^^^ as a trial function and adjust the parameters Coi^OiXo to cancel 
the contribution of the constant terms. Then we find that 

ip = iPo+i{G} = ^0 + ^ {e^' - eo^') = ^e". m 
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Actually, we can keep the constant term in the function T and find 

T = T^+ X{g2} = - ^ InC/Co = ^1 - ^ InC, (46) 

with the new parameter T\ = J"o + ^ In^o- 

Next, we sec that, fortunately, <1>('(/;^/,^) = <I>(Cq/4) is constant. Therefore we may continue 
applying the MIE procedure using an arbitrary DSG potential 

- e^^^{Cll\) / cg/4+./2_,-Cg/4+./2>._ 

-Co74 + i./2^ • ^ 

(48) 

The final integration in the first MIE iteration provides the improved function Q that can be used 
in the next iteration: 

_ e^^Co$'(CoV4) _^_cg/4+./2-i/2_ (49^ 



-Co74 + - 1/2 

(50) 

Of course, to avoid singularities, we suppose that i^7^Co/2 + l, v ^ Cq/2. 

We see that the new function Q has the same form g{C) = CiC~V2^ where Ci is a new constant, 

^ Co$-(CoV4)[2^-Cg] 
2$(CoV4)[2(^.-l)-C2]- 

If we choose Co satisfying the equation Ci = Co, the above expressions will give an exact solution. 
It can be useful for the investigation of the trajectories in vicinity of some singular points. 



4 Conclusion and outlook 

In a set of previous works on the topic of affine generalization of gravity we discussed the general 
properties of the theory and derived a very useful DSG interpretation. Here we tried to advance 
toward investigation of the dynamics provided by the theory. Thus the specific properties of the 
DSG dynamical system were revealed. Compared to the closely related but different DGS models, 
the DSG configuration appeared to be one of the most favorable for the analysis, possessing a 
linear inner symmetry and being well suited for the iterative MIE procedure. 

We suppose that this is not an accident but comes from the fundamental geometric origin of 
the affine gravity, the same as of the DGS models obtained by various dimensional reductions. 
In fact, such relations are usually hidden and can hardly be explicitly demonstrated on the level 
of Lagrangians. So, for the next step we will try to choose a somewhat different approach that 
would allow to establish a closer relation between the global (topological) properties of the set of 
all solutions and the inner symmetries of the DGS and DSG models. 

In conclusion, we wish to emphasize that the methods discussed above can be applied outside 
the context of the affine generalization of gravity and can provide, as we hope, useful tools for a 
wide class of models with scalar fields. 
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APPENDIX A: Dynamical systems with linear symmetries 

Consider a vanishing Hamiltonian for n dynamical variables g*, which is bilinear in momenta: 

H = A'^{q)piPj + U{q)=0. (A-1) 

Let us introduce in the coordinate space an affine transformation that defines the vector field 

v' = B'jq^ + c\ where Sj, = const. (A-2) 

For simplicity we will denote the vectors like {q^} by the bold font, q. The Euclidean scalar 
product, pv = (p, v) = PiV^, is a so-called Hamiltonian of the vector field v. Its Poisson bracket 
with the original Hamiltonian H obviously is 

{H, pv} = 2(p, BAp) - d,H, (A-3) 

where dy = v^d/dq^ is a derivative along the vector field v. 

Consider the linear vector fields the matrices of which satisfy the conditions: 

(p,B^p) = ^(p,Ap) ^ {B,A} = BA + AB = \A + B, (A-4) 

where B is an arbitrary antisymmetric matrix. Since H = (p, j4p) + [/ = 0, the above relation 
makes it possible to express the Poisson bracket (|A-3p as 

{H, pv} = 2(p, BAp) - dyH = A(p, Ap) - (p, d.Ap) - d^U = -\U - (p, d^Ap) - d^U. (A-5) 

It follows that the Poisson bracket of the two Hamiltonians, H and pv, vanishes when 

{B,A{q)}=XA{q)+B{q), dyA'^{q)=0, d,U{q) = -XU{q). (A-6) 

So, in case of the corresponding symmetries of the kinetic matrix and the potential term, dynamical 
systems (jA-ip admit the first integrals, which are linear in momenta and coordinates and are the 
Hamiltonians of the linear vector fields v = Bc[ + c. Their algebra reads 

{pv, pv'} = pv", where v" = [B' , B]q + {B'c - Be'), (A-7) 

where the square brackets denote the matrix commutator. 

One can easily resolve the conditions ()A-6p for the subgroup of transformations, v = Aq + c, 
which represent dilatations and translations: 

^^^■(q)=1>,,(q^), (A-8) 
A = : C/(q) = $i(q±), q± = jc^q - (q, c)c; (A-9) 
A/0: C/(q) = |vr2$2(v/|v|). (A-10) 

Here are the arbitrary functions of their arguments, which are now not independent: q^ 

is orthogonal to c, and v/|v| has unit length. 

APPENDIX B: Linear integrals in DGS 

Now let us return to the system ([2]) and try to construct integrals that are linear in momenta. For 
simplicity, we consider only the symmetries corresponding to the constant vector fields, v = c. 
Then the first relation in the conditions ()A-6p is automatically fulfilled for A = 0, since B is now 
zero matrix. 
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The remaining conditions can be satisfied in the following way. Consider first the dependence 
of the kinetic matrix only on the dilaton, Aij = ((/?), where the /i-dependence was removed by 
using the new variable F = h/h. The coordinate vector is now q = (F, ^i, . . . ,iPn)- Thus a 
possible ansatz has the form c = (a, 0, . . . , fi^), where a, /3* are constants. To compactify the 
lengthy formulas let {/3*} and {V^*} be treated as vectors /3, ^ in M^. 

All such constant fields, c, obviously satisfy dcAjk{q) = ddAjk{^p)/dq^ = 0. The solution for 
the potential term for such vector field is already obtained in (|A-9p . Construct first the vector q_L 
— the part of q that is orthogonal to c: 

= |/3|2f - a^^f, q^ = (a^ + qi = (a^ + - {aF + i = 1..N . (B-1) 

The potential can be an arbitrary function of q_L components: ^{F; {p;ipi, . . . , ij^^) = ^{c{±). Yet, 
in the actual gravitational potential, the metric enters as . So we can take the component 
as the argument of the exponent function. The metric can be excluded from the remaining 

components q\ if we use the scalar products q^ji'^^^ = (a^ + |/3p)/3jj™'''0 with A^ — 1 linearly 

independent vectors orthogonal to /3: /3/3|l"^ = 0, m = — 1. Then we have the following 
expression for the potential: 

e^Ui^-m = e^-^^^^"/l^^l'A(v9;{/3f VI), m = (B-2) 

where X is an arbitrary function of A^ arguments. The corresponding linear integral is 

pc = 00 + 2Zij {ip)/3'ip^ = P . (B-3) 

In the same way we can add the dependence of Zij on -i/^j^ , k = 1, .., AT. Then we should consider 
the constant vector /3 with vanishing those i^-t^ coordinates, actually belonging to M^~^. And 
the potential satisfying our conditions will be 

e^U{^-i;) = e^-'^^''/\^\" X{^-{'iPi^],{0'^^^]), k = l..K, m = l..N-K-l, (B-4) 

where now ^^^^ are A^ — A' — 1 arbitrary linearly independent constant vectors, belonging to the 
the orthogonal complement to /3 in M^~^ (they also have vanishing i^-th components). The first 
integral will be 

pc = a0 + 2Zi,-(c^;{^,J)/3*^J' = P. (B-5) 
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